# simulate a dataset with continuous data
y <- rdag(1000, 15, 0.3)
tru <- y$G
x <- y$x
mod <- pc.con(x)
mod$runtime
b <- pc.or(mod)
plotnetwork(tru)
dev.new()
plotnetwork(b$G)
dev.off()
plotnetwork( dag2eg(tru) ) ## essential graph
dev.new()
plotnetwork(b$G)
Run the code above in your browser using DataLab